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A Preliminary  Investigation  Into  Using  the  Sun 
as  a Source  for  G/T  Measurements 

Wi 1 1 i am  C . Daywitt 

Electromagnetic  Fields  Division 
National  Bureau  of  Standards 
Boulder,  Colorado  80303 

This  report  describes  a preliminary  investigation  into  determining  the  solar 
flux  density,  the  atmospheric  correction  factor,  and  the  star  shape  correction 
factor  for  use  in  G/T  measurements  above  5 GHz.  An  estimate  of  errors  is  also 

included.  Preliminary  results  show:  an  improved  algorithm  for  determining 

diffusive  and  refractive  attenuation;  a viable  technique  for  estimating  the  solar 
flux  density  from  daily  AFGL  flux  density  measurements  and  a centimeter/millimeter 
wave  spectrum  function;  and  the  possibility  of  reducing  star  shape  correction  factor 
errors  by  use  of  an  effective  solar  rf  diameter. 

Key  words:  atmospheric  correction  factor;  earth  terminal  measurement  system;  error 
analysis-,  G/T;  solar  flux  density;  star  shape  correction  factor. 


1.  Background  and  Introduction 

The  Electromagnetic  Fields  Division  of  the  National  Bureau  of  Standards  has  successfully  made 
precision  G/T  (gain  to  system  noise  temperature)  measurements  [1]  on  low  noise,  large  diameter 
(11  m to  28  m)  antenna  systems  using  the  radio  star  Cassiopeia  A,  and  on  medium  size  (4  m to  6 m) 
antenna  systems  using  the  moon.  Although  the  moon  is  an  attractive  source  to  use  with  these 
smaller  systems  because  of  its  predictable  radio  output  [2],  in  many  cases  [3]  the  system  noise  is 
so  much  greater  than  the  received  lunar  output  that  the  measurement  sensitivity  is  too  low  to  be 
useful.  Since  the  solar  angular  diameter  is  about  that  of  the  moon  and  its  output  is  considerably 
greater,  the  sun  provides  an  answer  to  the  sensitivity  problem.  Unfortunately,  its  flux  density 

and  brightness  distribution  are  variable,  and  even  with  this  preliminary  study  a number  of 
investigations  are  still  needed  to  determine  the  correct  G/T  measurement  error.  This  report 
describes  an  effort  to  determine  the  feasibility  of  using  the  daily  Air  Force  Geophysics 
Laboratory  (AFGL)  flux  density  values  [4]  measured  at  the  Air  Force  Sagamore  Hill  Observatory  in 
Hamilton,  Massachusetts.  The  results  pertain  to  frequencies  above  5 GHz  for  reasons  that  will  be 
made  clear  later. 

The  antenna  system  G/T  can  be  expressed  as 

G _ 8-nk  ( Y - 1) 

T 2 

x Sk1k2 

where  k is  Boltzmann's  constant,  Y is  the  ratio  of  the  power  measured  with  the  sun  in  the  - a 

mainbeam  to  the  power  measured  with  the  antenna  pointed  at  the  cold  sky,  \ is  the  v.  • ••  • 

wavelength,  S is  the  solar  flux  density  incident  on  the  earth's  atmosphere,  and  k|  a"  1 • , 
atmospheric  transmission  and  star  shape  correction  factors  respectively  [5]. 
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The  atmospheric  transmission  factor  accounts  for  the  reduction  in  the  received  flux  and  is 
due  to  three  atmospheric  phenomena  [6],  gaseous  attenuation,  refractive  spreading  of  the  incident 
wavefront,  and  incoherent  tropospheric  scattering.  It  does  not  include  attenuation  and  scattering 
from  clouds,  rain,  fog,  hail,  or  snow  since  the  precision  measurements  envisioned  here  are  assumed 
to  be  performed  in  clear,  stable  weather.  If  A is  the  total  zenith  attenuation  in  decibels  due  to 
these  three  phenomena,  then 


kj  . i0-'vl0s1,w  (2) 

where  0 is  the  elevation  angle  to  the  source  relative  to  the  antenna  location  on  the  earth's 
surface.  Details  of  a simplified  calculation  for  the  gaseous  component  of  A are  described 
elsewhere  [7], 

An  algorithm  derived  in  an  earlier  analysis  [5]  for  the  zenith  attenuation  due  to  the 
tropospheric  scatter  of  the  energy  from  a point  source  radio  star  and  intended  for  the  1 GHz  to  10 
GHz  frequency  range  diverges  above  10  GHz.  Since  that  analysis,  further  measurement  data  [8]  has 
come  to  light  which  provides  a more  detailed  picture  of  this  attenuation.  The  resulting  curves 
derived  from  the  data  and  which  include  the  refractive  attenuation  are  shown  in  figure  1.  The 
first  set  of  data  [6]  were  taken  at  4.17  GHz  and  6.39  GHz  with  a 22  m antenna,  and  the  second  set 
[8]  at  15.5  GHz  and  31.6  GHz  with  a 7 m antenna.  The  dashed  curves  through  the  two  sets  of  data 
points  are  parallel,  and  suggest  a family  of  straight  lines  for  other  antenna  diameters.  A linear 
interpolation  in  antenna  diameter  between  the  two  dashed  curves  produced  the  linear  portions  of 
the  curves  shown.  Since  theory  [6]  indicates  that  the  attenuation  decreases  with  decreasing 
frequency,  the  curves  were  extended  accordingly.  When  the  angular  diameter  of  the  source  is 

finite,  the  attenuation  is  reduced  [6],  decreasing  as  the  ratio  of  source  diameter  to  antenna 
beamwidth  increases.  The  equation  in  figure  1 descibes  this  effect  [5]  where  d is  the  source 
diameter,  is  the  half-power  beamwidth,  AQ  is  the  zenith  attenuation  for  a point  source  given  by 
the  curves  in  the  figure,  and  A is  the  corresponding  attenuation  for  the  finite  source  diameter. 
An  earlier  investigation  [5]  indicates  that  the  total  error  in  these  attenuation  curves  and  the 
non-point  source  correction  is  probably  in  the  neighborhood  of  50%. 

The  star  shape  correction  factor  accounts  for  variation  of  the  solar  brightness  distribution 
across  the  antenna  beam  pattern,  and  vanishes  as  the  ratio  of  the  angular  source  diameter  to  the 
antenna  half-power  beamwidth  vanishes.  In  terms  of  a distribution  B and  a normalized  power 
pattern  Pn, 


/ B(q)P  (r)  do 

k2  = / B(r)  dR  (3) 

where  R represents  the  two  antenna  pointing  angles  and  dR  is  the  differential  solid  angle.  The 
integral  in  the  denominator  is  equivalent  to  the  flux  density  S. 

Sections  2 and  3 present  the  results  of  a preliminary  investigation  into  the  determination  of 
the  flux  density  S and  the  star  shape  correction  factor 
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2.  Determination  of  the  Solar  Flux  Density  Above  5 GHz 
2.1  Strategy 


Solar  radio  emission  is  the  result  of  three  superimposed  sources  of  radiation  [9],  the  quiet, 
slow,  and  burst  components.  The  quiet  sun  emission  comes  from  unlocalized,  thermal  sources  in  the 
solar  atmosphere  (chromosphere  and  corona),  and  is  constant  over  periods  of  months  to  years.  The 
slow  component  also  has  a thermal  origin  and  originates  in  local  sources  in  the  vicinity  of 
sunspots  and  chromospheric  plagues.  These  partially  polarized  sources  have  periods  of  days, 
weeks,  or  months.  The  burst  components  are  the  results  of  thermal,  synchrotron,  and  plasma 
radiation  from  transient  sources  such  as  flares,  and  are  partially  polarized.  Their  recurrence 
frequency  and  period  are  small  for  the  quiescent  portion  of  the  solar  cycle,  increasing  to  two  or 
more  per  a 10-hour  viewing  day  during  the  active  portion,  which  suggests  the  possibility  of  making 
G/T  measurements  using  the  sun  without  encountering  a burst  during  the  process.  The  fact  that  the 
quiet  and  slow  components  have  periods  of  greater  than  a day  further  suggests  the  possibility  of 
using  an  average  of  two  successive  daily  AFGL  flux  density  values  to  determine  the  flux  at  the 
time  of  the  G/T  measurement  taking  place  between  the  two  AFGL  measurements. 

Figure  2 depicts  the  average  solar  flux  density  spectrum  incident  on  the  earth's  atmosphere 
[10]  with  an  insert  [11]  indicating  the  duration  of  the  burst  spectrum  from  10  MHz  to  30  GHz.  The 
quiet  component  of  the  solar  emission  dominates  above  5 GHz,  at  which  frequency  most  of  the  burst 
components  have  disappeared  and  the  slow  component  is  falling  off.  The  insert  shows  that  the 
duration  of  the  remaining  burst  continuum  is  short.  Thus,  the  solar  emission  above  5 GHz  is 
relatively  stable  from  day  to  day  as  long  as  the  type  IV  bursts  are  avoided,  again  suggesting  a 
G/T  measurement  scheme  using  AFGL  flux  density  values. 

Superimposed  on  the  graph  of  figure  2 are  three  sets  of  hashmarks  representing  monthly 
averages  of  AFGL  flux  density  measurements  (with  atmospheric  attenuation  removed)  [4]  from  January 
1979  to  May  1982.  The  measurement  values  approach  the  average  quiet  sun  value  between  the 
"sunspot  (s.s)  max"  and  the  "s.s.  min"  curves  as  the  frequency  increases.  Thus,  above  the  highest 
AFGL  measurement  frequency  (15.4  GHz),  apparently  the  average  of  the  quiet  sun  curves  can  be  used 
to  determine  the  flux  density  incident  on  the  earth's  atmosphere.  Following  is  an  approximate 
equation  giving  this  average  that  was  derived  by  least-square  fitting  the  average  data 
corresponding  to  the  quiet  sun  curves  in  the  figure  and  in  table  1 (from  reference  [10]): 


The  flux  density  S and  the  frequency  f are  the  same  units  as  the  figure.  The  following  equations 
can  be  used  to  extrapolate  between  the  AFGL  measurement  frequencies  [12]: 


log  S = 1.20  + 1.10  (log  f)  + 0.179  (log  f)2. 


(4) 


where 


1 n Cf /f ! ) 


(6) 


S has  the  same  units  (W/M2/Hz)  as  the  graph,  and  f i and  f£  correspond  to  Sj 
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and  S^. 


respectively.  The  subscripts  denote  the  two  AFGL  values  of  flux  density  and  frequency  between 
which  the  flux  density  S at  frequency  f is  desired. 

The  preceding  discussion  suggests  the  following  strategy  for  determining  the  solar  flux 
density  above  5 GHz:  1)  if  the  G/T  measurement  frequency  is  4.995  GHz,  8.8  GHz,  or  15.4  GHz,  use 
an  average  of  two  adjacent  daily  AFGL  flux  values;  2)  if  the  measurement  frequency  is  between  two 
AFGL  frequencies,  use  the  average  determined  in  1)  for  both  AFGL  frequencies  and  translate  the 
result  to  the  desired  frequency  via  eqs  (5)  and  (6);  and  3)  if  the  frequency  is  above  15.4  GHz, 
use  eq  (4)  to  determine  the  flux  density.  Since  the  reported  AFGL  flux  density  values  are 

corrected  for  bursts,  it  is  only  then  necessary  to  insure  that  the  G/T  measurement  did  not  take 
place  duri  ng  a burst . 

To  illustrate,  suppose  a series  10  GHz  G/T  measurements  were  performed  in  Boulder,  Colorado, 
at  various  times  between  7 a.m.  (1400  u.t.  (universal  time))  and  5 p.m.  (2400  u.t.  local  standard 
time  on  October  17  , 18,  19,  1983.  The  appropriate  AFGL  flux  density  values  (in  units  of  10“^ 
W/m/Hz)  that  bracket  10  GHz  on  these  days  are  isolated  in  the  box  of  figure  3 for  the  month  in 
question  [4].  The  "C"  value  from  eq  (6)  for  the  frequencies  indicated  (8.8  GHz,  10  GHz,  15.4  GHz) 
turns  out  to  be  0.7716  (with  f^  = 8.8,  f = 10,  and  ^ ~ 15.4).  Using  this  value,  the  flux 

densities  in  the  box,  and  eq  (5),  the  10  GHz  flux  density  values  at  the  AFGL  measurement  time 

(approximately  1630  u.t.)  on  the  17th,  18th,  19th,  and  the  20th  are  334,  321,  330,  and  326, 

respectively.  Linear  interpolation  between  any  appropriate  two  of  these  values  will  then  yield 
the  correct  flux  value  for  a specific  Boulder  G/T  measurement  time.  Figure  4 is  a reproduction  of 
another  page  from  the  "Solar-Geophysical  Data"  report  that  indicates  the  time  and  size  (in  flux 
density  unito)  of  any  bursts  that  have  taken  place  during  the  days  in  which  the  G/T  measurements 
were  being  performed.  The  data  in  the  boxes  show  the  larger  microwave  bursts  taking  place  during 
the  G/T  measurement  times  in  question  (1400  u.t.  to  2400  u.t.),  and  thus  any  G/T  measurements 
performed  during  these  bursts  chould  be  discarded.  The  smaller  bursts  shown  in  the  figure, 
although  somewhat  arbitrarily  excluded,  have  a less  significant  effect  on  the  G/T  measurement. 

2.2  Flux  Density 

For  simp'Mcity  the  preceding  illustration  ignores  atmospheric  attenuation  and  the  true  earth- 
sun  distance  at  the  time  of  thp  G/T  measurements.  A complete  flux  density  determination  between 
AFGL  frequencies  proceeds  as  follows:  1)  calculate  the  flux  densities  S-^  and  S2  corresponding  to 
the  appropriate  two  AFGL  frequencies  for  the  G/T  measurement  time  between  the  two  AFGL  daily 
measurements;  2)  adjust  and  S2  for  the  true  earth-sun  distance  (the  values  reported  in  figure  3 

are  normalized  to  1 astronomical  unit);  3)  remove  the  effects  of  atmospheric  attenua'^on  from 
and  S2  via  eq  (2);  4)  use  eqs  (5)  and  (6)  to  obtain  the  flux  density  S incident  on  tlr  " ch's 
atmosphere  at  the  G/T  measurement  frequency;  and  5)  apply  eq  (2)  a second  time  at  the  measurement 
frequency  to  obtain  the  flux  density  at  the  location  on  the  earth's  surface  where  the  G/T 
measurement  is  performed.  Above  15.4  GHz  eq  (4)  can  be  used  to  estimate  the  flux  density  incident 
on  the  atmosphere,  after  which  step  5 is  applied. 

Prior  to,  but  not  afterward,  January  1980,  AFGL  removed  the  effects  of  atmospheric  attenua- 
tion from  their  measurements  and  reported  flux  density  values  incident  on  the  atmosphere. 
Therefore,  the  flux  density  values  reported  in  figure  3 contain  the  atmospheric  attenuation. 
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To  illustrate,  suppose  a G/T  measurement  is  performed  at  an  antenna  site  in  Boulder  on 

October  17,  1983,  as  the  sun  transits  the  local  meridian  (1848  u.t.  [13]).  The  AFGL  flux  density 
values  measured  earlier  on  the  17th  and  on  the  following  day  during  meridian  transit  (1630  u.t. 
[13])  at  the  Sagamore  Hill  Observatory  are  reported  (fig.  3)  as  282  and  268  at  8.8  GHz,  and  as  588 
and  587  at  15.4  GHz.  By  a linear  extrapolation  in  time  between  these  daily  measurements,  flux 
density  values  at  1848  u.t.  of  280.6  and  587.9  are  obtained  for  8.8  GHz  and  15.4  GHz, 

respectively.  The  true  earth-sun  distance  (in  units  of  1 astronomical  unit)  at  1848  u.t.  on  the 
17th  was  0.996  [13].  Thus,  the  true  flux  density  values  at  the  G/T  measurement  time  are  a factor 

of  1 .008  (0.996-^)  larger  than  the  values  reported  by  AFGL,  or  282.8  and  592  .6  for  8.8  GHz  and 

15.4  GHz,  respectively. 

The  average  solar  declination  during  the  AFGL  and  G/T  measurements  was  a negative  9.4  degrees 
[13].  Combined  with  the  Sagamore  Hill  latitude  of  40.6°  north,  this  leads  to  a solar  elevation 
angle  at  that  site  of  38.0°.  The  zenith  values  for  the  gaseous  component  of  the  atmospheric 

attenuation  in  October  at  the  observatory  are  [7]  0.042  dB  and  0.080  dB  for  8.8  GHz  and  15.4  GHz, 
respectively.  The  zenith  values  for  the  diffusive  and  refractive  attenuation  for  the  average 

solar  diameter  of  0.58°  [13]  during  the  measurements  are  (fig.  1)  0.012  dB  at  8.8  GHz  and  0.015  dB 
at  15.4  GHz.  Combining  these  with  the  gaseous  components  gives  a total  zenith  attenuation  of 
0.054  dB  at  8.8  GHz  and  0.095  dB  at  15.4  GHz.  These  latter  values  with  the  38.0°  declination  lead 
via  eq  (2)  to  atmospheric  correction  factors  of  0.980  at  8.8  GHz  and  0.965  at  15.4  GHz.  Dividing 
the  previous  flux  density  values  by  these  correction  factors  leads  to  288.6  and  614.1  at  8.8  GHz 
and  15.4  GHz,  respectively,  for  the  flux  densities  incident  on  the  atmosphere.  Equations  (5)  and 

(6)  can  now  be  used  to  obtain  an  estimate  of  the  10  GHz  flux  density  and  results  in  a value 

342.9.  If  the  Boulder  site  is  at  40°  north  latitude,  the  solar  elevation  angle  there  is  40.6°. 
Combining  a gaseous  attenuation  of  0.029  dB  [7]  with  a diffusive  and  refractive  attenuation  of 
0.011  dB  (fig.  1)  yields  a correction  factor  of  0.986  for  a 40.6°  elevation  angle,  resulting  in  a 

value  of  338.1  for  the  flux  density  at  the  Boulder  site. 

2.3  Errors 

The  total  error  in  determining  the  flux  density  from  the  AFGL  measurements  is  the  sum  of  four 

parts:  the  AFGL  measurement  error;  the  error  in  using  eqs  (5)  and  (6)  to  extrapolate  between  AFGL 

frequencies;  the  error  in  estimating  the  atmospheric  transmission  factor;  and  the  error  in 
estimating  the  flux  density  by  the  linear  time  interpolation  between  the  daily  AFGL 
measurements.  This  latter  error  for  frequencies  above  the  highest  AFGL  frequency  (15.4  GHz)  is  at 
present  taken  as  the  maximum  error  (halfway)  between  the  highest  AFGL  frequencies,  8.8  GHz  and 

15.4  GHz. 

The  AFGL  error  for  measuring  the  flux  density  is  5%  [12]. 

An  error  of  1%  is  assigned  to  the  uncertainty  in  extrapolating  between  AFGL  frequencies,  and 
was  determined  by  comparing  the  8.8  GHz  average  value  (average  of  the  AFGL  hashmarks  in  figure  2) 
with  the  8.8  GHz  value  obtained  via  eqs  (5)  and  (6)  and  the  average  AFGL  values  at  4.995  GHz  anc 

15.4  GHz. 

The  total  (gaseous  + diffusive  + refractive)  zenith  attenuations  for  Boulder  at  10  GHz  r : 
Sagamore  Hill  at  8.8  GHz  and  15.4  GHz  are  0.040  dB,  0.054  dB,  and  0.095  dB,  respectively, 
these  values  are  uncertain  by  50%  [5],  leading  to  transmission  coefficient  errors  of  0.7  , . 
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and  1.8%,  respectively.  The  overall  transmission  coefficient  error  is  a weighted  sum 


(0 .7+cxl .0+( 1-c )xl .8)  of  these  separate  errors  with  a value  of  approximately  2%  for  the  value  of  C 
determined  in  section  2.1. 

Figure  5 is  a plot  of  a number  of  hypothetical  smoothed  flux  density  curves  between  the  AF6L 
values  of  S(t^)  and  S(t2)  measured  during  meridian  passage  at  time  tj  and  time  t2  the  following 
day  (t^  - t^  + 24  hours).  The  linear  curve  whose  equation  is  shown  in  the  figure  was  previously 
used  to  estimate  the  flux  density  at  time  t between  the  AFGL  measurement  times  t^  and  t2,  where  h 
is  the  relative  forward  interval  between  t^  and  t.  S(t)  is  the  flux  density  at  t for  the  true 
(but  unknown)  curve,  and  as  is  the  estimated  standard  deviation  of  the  possible  values  at  t for 
the  probability  density  curve  shown.  The  forward  relative  interval  h is  defined  in  the  figure. 
AS  is  defined  to  be  the  difference  between  the  linear  approximation  and  the  unknown  flux  density 

S(t),  and  is  therefore  a random  variable  whose  estimated  standard  deviation  a is  a measure  of  the 

A 

error  in  using  the  linear  approximation.  A simple  analysis  of  the  situation  depicted  in  the 
figure  shows  that  the  two  standard  deviations  are  related  by  the  equation 


where  p,  is  the  estimated  autocorrelation  function  [14]  between  measured  flux  density  values  24 
■L  h 1 h 

hours  (=  t^  - t^)  apart.  and  p^-  are,  respectively,  approximations  to  the  autocorrelation 

functions  between  flux  density  values  t - tj  and  t2  - t apart.  The  flux-density  variations  are 
assumed  to  be  sufficiently  stationary  for  eq  (7)  to  give  some  indication  of  the  time  interpolation 
error . 

Figure  6 shows  a plot  of  the  8.8  GHz  AFGL  flux  density  values  measured  from  January  through 
June  1979.  The  change  of  levels  at  LI,  L2,  and  L3,  and  the  slopes  at  SI,  S2,  and  S3  indicate  that 
the  variation  is  not  stationary  and  should  possibly  be  modeled  [14]  as  an  ARIMA  (p ,1  ,q ) or  an 
ARIMA  (p,2,q)  process  to  estimate  the  standard  deviation  of  the  residual  day-to-day  fluctuations 
of  interest.  However,  for  a preliminary  (and  more  conservative)  estimate  of  os , the  data  shown 
was  assumed  to  be  stationary,  resulting  in  a value  of  24  for  os  and  a value  of  0.9  for  p-^.  With 
these  values  eq  (7)  is  approximated  by 


This  equation  shows  a maximum  of  5.5  halfway  between  tj  and  t2,  and  is  assumed  to  hold  for  the 
entire  range  from  5 GHz  to  15.4  GHz.  For  the  illustration  in  section  2.2,  h has  the  value  of 
0.096  ( ( 18 .8-16 .5 ) / 24 ) , yielding  an  estimated  standard  deviation  from  eq  (8)  of  3.2  or  1% 

(3.2/337),  or  a 3a  error  of  approximately  3%. 

Collecting  the  previous  errors  leads  to  a total  error  of  11%  (3a)  in  determining  the  solar 
flux  density  for  the  illustration  presented  in  section  2.2.  The  error  at  lower  elevation  angles 
and/or  higher  frequencies  will  be  greater,  primarily  because  of  increased  gaseous  attenuation. 

3.  Star  Shape  Correction  Factor 

For  antenna  half-power  beamwidths  (HPBW)  greater  than  about  five  angular  source  diameters, 
examination  of  eq  (9)  shows  that  ratio  in  eq  (3)  can  be  replaced  by  unity  with  an  error  less  than 


- 2a^[l  - h + h^  + h(l-h)p^  - (l-h)p^  - hpj  ^] 


(7) 


(8) 
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l‘t. . To  use  the  cel estial -source  method  of  measuring  G/T  accurately  below  five  source  diameters, 
however,  an  estimate  of  the  variation  of  the  source  brightness  distribution  across  the  antenna 
mainbeam  indicated  in  the  numerator  of  eq  (3)  must  be  made.  For  sources  with  stable  and 

predictable  distributions  and  Gaussian  mainbeams,  this  can  be  accomplished  by  replacing  the  ratio 
in  eq  (3)  by  a simple  formula  depending  only  upon  the  ratio  of  HPBW  to  source  diameter. 
Furthermore,  if  the  source  diameter  is  properly  chosen,  the  formula  will  provide  an  accurate 
estimate  of  the  correction  factor  for  any  HPBW/di ameter  ratio  above  unity  [15,2].  As  will  be  seen 
there  is  reason  to  expect  that  this  technique  will  also  work  for  the  sun,  even  though  the  solar 
brightness  distribution  is  continuously  changing.  Errors  in  the  star  shape  correction  factor  have 
not  yet  been  investigated. 

For  a source  resembling  a circular  disk  and  a Gaussian  mainbeam,  eq  (3)  leads  to 


where 


2 


x 


2 

1 n2  . 


(9) 


(10) 


d is  the  effective  rf  angular  diameter  of  the  source,  and  is  the  antenna  HPBW.  The  diameter  is 
chosen  to  minimize  the  difference 


1 - e“x  / B (q ) Pn(Q)  dQ 
x2  / B(q)  dQ 


(ID 


for  eH/d  greater  than  unity.  The  same  expressions  (eqs  (9)  and  (10))  hold  for  an  elliptical  disk 

O 

with  dc  replaced  by  the  product  of  the  major  and  minor  source  diameters.  No  significant  reduction 
in  the  difference  in  eq  (10)  is  achieved  with  this  latter  model,  however,  even  though  the  solar 
brightness  distribution  shows  a decided  ellipticity  (fig.  7,  [9]). 

Figure  2 shows  that,  except  for  the  burst  components  which  are  assumed  to  be  missing  during 
the  G/T  measurement,  the  quiet  sun  emission  plays  an  increasingly  dominant  role  in  the  total 
emission  above  5 GHz  where  the  slow  component  begins  to  fall  off.  Furthermore,  observations  of 
the  solar  brightness  distribution  in  the  decimeter,  centimeter,  and  millimeter  wave  frequency 
regions  reveal  that,  although  the  effective  brightness  temperature  varies  with  the  solar  cycle, 
the  distributions  remain  relatively  constant  [9].  It  is  this  constant  part  of  the  brightness 

distribution  that  will  presumably  be  amenable  to  the  technique  of  the  previous  paragraph,  wit’! 
fluctuations  in  the  result  coming  from  the  more  variable  slow  emission  component.  The  degree  of 
success  will  depend  upon  how  minor  these  fluctuations  turn  out  to  be,  but  hopefully  at  least  sir- 
accuracy  in  k£  will  be  gained  below  five  source  diameters.  This  question  can  only  be  de-rv: 
after  a number  of  distributions  at  various  frequencies  and  at  various  times  during  the  solar  cy<  1- 
have  been  convoluted  with  the  antenna  beam  pattern.  To  date  this  has  not  been  done. 

Since  the  effective  diameter  or  diameters  for  the  solar  distributions  have  not  yet 
determined,  the  question  of  what  diameter  d to  use  in  eq  (10)  remains.  A glance  at  f 1 , r - • 
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that  the  polar  profiles  tend  to  average  out  as  a constant  up  to  the  optical  diameter,  with  the 
equatorial  profile  averages  extending  somewhat  beyond.  This  situation  implies  that  the  effective 
rf  diameter  will  extend  beyond  the  optical  diameter,  which  in  turn  will  depend  upon  the  results 
from  convolving  the  brightness  distributions  with  the  Gaussian  beam  pattern.  Based  on  a circular 
disk  model  of  the  quiet-sun  atmosphere  with  a constant  chromospheric  temperature  of  3 x 10^ 
kelvins  and  a coronal  temperature  of  10®  kelvins,  Smerd  [16]  calculated  a brightness  temperature 
distribution  and  compared  it  to  the  apparent  brightness  temperature  derived  from  the  same  model. 
The  apparent  temperature  is  defined  as  the  temperature  of  a blackbody  having  the  same  diameter  as 
the  optical  disk  and  radiating  the  same  apparent  flux  density.  The  results  are  shown  in  figure  8 
where  the  ratio  of  two  particular  brightness  temperature  contours  to  the  apparent  temperature  are 
plotted.  The  upper  curve  is  where  the  contour  temperature  falls  to  IX  of  the  apparent  temperature 
and  the  lower  curve  is  where  it  falls  to  10%.  As  expected  the  ratios  approach  unity  with 
increasing  frequency  as  the  origin  of  the  radiation  is  increasingly  confined  to  the 
chromosphere.  The  dashed  line  approximating  the  10%  curve  and  ending  at  30.3  GHz  will  be  used 
temporarily  to  define  the  effective  diameter  for  eq  (10).  Thus, 


d = d 

o 


1 .240-0.162  log  f,  5<f<30.3 
1 , f>  30.3 


(12) 


where  dQ  is  the  optical  diameter  obtained  from  the  Astronomical  Almanac  [13]  at  the  time  of  the 
G/T  measurement.  For  the  illustration  in  section  2.2,  the  Boulder  antenna  diameter  is  4.3  m (fig. 
1)  yielding  a 10  GHz  HPBW  of  0.49°.  At  the  measurement  time  the  optical  solar  diameter  was  0.54° 
[13]  which  yields  an  rf  diameter  via  eq  (12)  of  0.58°.  Equations  (9)  and  (10)  then  give  a star 
shape  correction  factor  of  0.64. 


4.  Discussion  and  Conclusions 

The  curves  in  figure  1 used  for  determining  the  diffusive  and  refractive  attenuation  fc-'  a 
point  source,  although  quite  rough  and  incomplete,  offer  an  improvement  over  the  previous 
algorithm  [5]  that  diverges  above  10  GHz. 

The  total  error  in  estimating  the  solar  flux  density  from  5 GHz  to  15.4  GHz  in  clear,  stable, 
weather  is  11%,  and  is  the  sum  of  5%  for  the  AFGL  measurements,  1%  for  extrapolati  ng  between  two 
AFGL  frequencies,  3%  for  estimating  the  flux  density  between  the  daily  AFGL  measurement  times, 
1.3%  for  removing  the  atmospheric  attenuation  from  the  AFLG  measurements,  and  0.7%  for  estimating 
the  atmospheric  attenuation  at  the  G/T  measurement  frequency.  Discarding  the  0.7%  error,  the  sum 
comes  to  about  8%,  and  represents  the  error  in  determining  the  flux  density  incident  on  the 

earth's  atmosphere.  Since  the  average  of  the  AFGL  hashmarks  in  figure  2 approaches  the  area 

between  the  "s  .s . max"  and  "s.s.  min"  curves  as  the  frequency  increases,  and  since  their  spread 
and  that  between  the  curves  decreases  with  frequency,  it  appears  that  an  8%  error  in  using  eq  (4) 
to  determine  the  solar  spectrum  above  15.4  GHz  incident  on  the  earth's  atmosphere  is  reasonable 
for  the  spectrum  shown  in  the  figure. 

The  brightness  distribution  of  the  quiet  sun  above  5 GHz  is  relatively  insensitive  to  the 

solar  cycle.  This,  coupled  with  the  fact  that  the  quiet-sun  component  also  dominates  above  this 
frequency,  implies  the  possibility  of  finding  an  effective  rf  diameter  (as  a function  of 

frequency)  that  minimizes  the  error  in  using  the  circular  disk  approximation  to  the  star  shape 
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correction  factor  in  eq  (3),  improving  G/T  measurement  accuracy  for  HPBW/source-di ameter  ratios 
between  one  and  five. 

Although  the  preceding  conclusions  appear  to  be  reasonable,  there  still  remains  a 
considerable  amount  of  work  to  be  completed  before  final  acceptance. 


The  majority  of  the  work  described  in  this  report  was  carried  out  in  1979  and  was  part  of  a 
larger  effort  to  design  and  construct  an  earth  terminal  measurement  system  capable  of  precision 
antenna  system  measurements.  Therefore,  much  of  the  direction  and  inspiration  for  this  work  came 
from  other  individuals  associated  with  the  project,  principally  David  Wait  and  John  Wakefield. 
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Figure  1.  Diffusive  and  refractive  attenuation  AQ  for  a point  source.  The  equation 
appearing  in  the  upper  left  corner  of  the  figure  is  used  to  estimate  the  attenuation  A 
for  a non-point  source  (d  * 0)  once  An  has  been  determined  from  the  curves. 


Solar  Flux  (W/m2/Hz) 


Figure  2.  Solar  radio 
reference  [10].) 


spectrum.  (The  symbols  are  explained  more  fully  in 
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The  observed  and  the  adjusted  Ottawa  fluxes  tabulated  above  are  the  "Series  C"  dally  values  reported  by 
the  Algonquin  Radio  Observatory,  Ottawa,  Ontario,  Canada.  The  letter  "A"  following  an  enfry  designates 
an  Interpolated  flux.  Numbers  in  parentheses  In  the  column  headings  denote  frequencies  In  MHz. 

Equipment  problems  produced  the  gaps  shown  here  In  the  Air  Weather  Service's  Sagamore  Hill  (SGMR)  obse> — 
vat i ons . 


The  International  and  American  sunspot  numbers  shown  above  are  preliminary  values. 


Figure  3.  Smoothed  solar  flux  density  (in  units  of  W/m^/Hz)  measured 
the  AFGL  Sagamore  Hill  observatory  (SGMR)  at  approximately  solar 
passage  each  day  and  adjusted  to  one  astronomical  unit.  The  box 
the  flux  densities  of  interest. 
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Reports  are  received  routinely  from  the  following  observatories: 


ATHN  = Athens 

HUAN  = 

Huancayo 

NAGO  = 

Nagoya 

°0TS 

= Potsdam 

BERN  = Berne 

IRKU  = 

1 rkutsk 

N06E  = 

Nobeyama 

SA0P 

= Sao  Paj 1 o 

BORD  = Bordeau 

IZMI  = 

IZMIRAN 

ONDR  = 

Ondrej  ov 

SGMR 

= Sagamore  Hill 

CR IM  = Cr imea 

KISV  = 

Ki slovodsk 

OTTA  = 

Ottawa 

TORN 

= Torun 

DW 1 N = Dw  i nge 1 ' 

DO 

KRAK  = 

Krakow 

PALE  = 

Pa  1 ehua 

TYKW 

= Toyokawa 

GORK  = Gorky 

LEAR  = 

Learmonth 

PEKG  = 

Pek i ng 

TRST 

= Trieste 

H IRA  = H I rai  so 

MANI  = 

Man i 1 a 

PENT  = 

Penticton 

UP  1C 

= Up  ice 

V0R0 

= Vorosh i 1 ov 

Explanation  of 

Type  Code: 

1 Simple  1 
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R i se 

30 

Post  Burst  Increase  A 

43  Onset  of  Noi 

se  Storm 

2 Simple  IF 
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Spike  25 
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31 

Post  Burst  Decrease 

44  Noise  Storm 
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3 Simple  2 

20 

Simple  3 26 

F a 1 1 

33 

Absorption 

45  Complex 
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Simple  3A  27 

R 1 se  and  Fall 
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46  Comp  1 ex_  F 
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Post  Burst  Increase  42 

Series  of  Bursts 

48  Major 
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Fal  1 F 
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21A  Simple  3A  GRF  240F  Rise  only  F 260  Fall  Only  31A  Post  Burst  Decrease  A 

2A  Simple  1AF  24P  Post  Rise  26F  Fall  F 32A  Absorption  A 

46F  Complex  F 

Remarks : 

QL  = Quality  ( 1 =poor  to  6=excel lent) 

ST  = Status  (1=real  time;  2=  final ; 3=correct ion ; 4=deletion) 

TYP=  Type  (1=noise  storm;2=r!se  In  base  level ; 3=minor;4=group;5=major;6=major  plus;7=Castel I i U-type  burst) 
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Figure  5.  Hypothetical  flux  density  curves. 


Figure  6.  Smoothed  solar  flux  density  at  8.8  GHz  measured  by  AFGl  f' 
January  to  August  1979. 
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Figure  7.  Quiet -sun  brightness  temperature  profiles  (in  units  of  10^  kelvins) 
at  various  frequencies.  The  solid  curves  are  equatorial  profiles,  and  the 
dashed  curves  are  polar  profiles.  The  values  of  e represent  the 
ellipticity  at  the  half-power  levels  indicated.  The  abscissa  is  in  units 
of  the  optical  solar  radius.  (After  Kundu  [9].) 
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Figure  8.  Effective  solar  rf  diameter  (d)  in  units  of  the  optical  diameter 
(dQ)  versus  frequency  for  two  ratios  (1%  and  10%)  of  contour  temperature 
to  apparent  temperature.  (After  Smerd  [16].) 


p 

Table  1.  Average  flux  density  S(W/m  /Hz)  versus  frequency  for  the  quiet  sun  radiation  in 
terms  of  the  base-10  logarithm  of  S.  (From  Allen  [10].) 


QUIET  SUN  RADIATION 

FREQUENCY  (GHz) 

Log  S 

10 

2.48 

20 

2.93 

50 

3.57 

100 

4.125 

200 

4.69 

500 

5.47 

18 


NBS-I  14A  REV.  2-g 


U.S.  DEPT.  OF  COMM. 

BIBLIOGRAPHIC  DATA 
SHEET  (See  instructions) 


1.  PUBLICATION  OR 
REPORT  NO. 

NBSIR  84-3015 


2.  Performing  Organ.  Report  No, 


3.  Publication  Date 

August  1984 


4.  TITLE  AND  SUBTITLE 

A Preliminary  Investigation  Into  Using  the  Sun  as  a Source  for  G/T  Measurements 


5.  AUTHOR(S) 


William  C.  Daywitt 


6.  PERFORMING  ORGANIZATION  (If  joint  or  other  than  NBS.  see  instructions) 

NATIONAL  bureau  of  standards 
department  of  commerce 

WASHINGTON,  D.C.  20234 


7.  Contract/Grant  No. 


8.  Type  of  Report  & Period  Covered 


9.  SPONSORING  ORGANIZATION  NAME  AND  COMPLETE  ADDRESS  (Street.  City.  State.  ZIP) 

Same  as  6 


10.  SUPPLEMENTARY  NOTES 


1 Document  describes  a computer  program;  SF-185,  FIPS  Software  Summary,  is  attached. 

11.  ABSTRACT  (A  200-word  or  less  factual  summary  of  most  significant  information.  If  document  includes  a significant 
bi  bl  iography  or  literature  survey,  mention  it  here) 

This  report  describes  a preliminary  investigation  into  determining  the  solar 
flux  density,  the  atmospheric  correction  factor,  and  the  star  shape  correction  factor 
for  use  in  G/T  measurements  above  5 GHz.  An  estimate  of  errors  is  also  included.  Pre- 
liminary results  show:  an  improved  algorithm  for  determining  diffusive  and  refractive 

attenuation;  a viable  technique  for  estimating  the  solar  flux  density  from  daily 
AFGL  flux  density  measurements  and  a centimeter/millimeter  wave  spectrum  function;  an 
the  possibility  of  reducing  star  shape  correction  factor  errors  by  use  of»an  effective 
solar  rf  diameter. 


12.  KEY  WORDS  (Six  to  twelve  entries;  alphabetical  order;  capitalize  only  proper  names;  and  separate  key  words  by  semicolon  sj 

atmospheric  correction  factor;  earth  terminal  measurement  system;  error  analysis; 
G/T;  solar  flux  density;  star  shape  correction  factor 


13.  AVAILABILITY 
| xl  Uni  imited 

| | For  Official  Distribution.  Do  Not  Release  to  NTIS 

[2ZI  Order  From  Superintendent  of  Documents,  U.S.  Government  Printing  Office,  Washington,  D.C. 
20402. 

|_Xj  Order  From  National  Technical  Information  Service  (NTIS),  Springfield,  VA.  22161 


14.  NO.  OF 

PRINTED  PAGES 

24 


15.  Price 

$7.00 


n 


h 


■ 


■ 


